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SUMMARY OF INFORMATION ON 
LOW-SPEED LATERAL-DIRECTIONAL DERIVATIVES DUE TO 
RATE OF CHANGE OF SIDESLIP /§ 

Paul L. Coe, Jr., A. Bruce Graham,* and Joseph R. Chambers 
Langley Research Center 

SUMMARY 

The present paper summarizes the experience obtained by NACA and NASA con- 
cerning the low-speed aerodynamic stability derivatives due to the rate of change of side- 
slip 0. The report also includes a comprehensive bibliography on this and related 
subject matter. 

The data show that large values of the acceleration-in-sideslip derivatives C n ^ 
and C are obtained for swept- and delta-wing configurations at high angles of attack. 
The physical flow phenomenon responsible for the 0 derivatives is associated with 
(1) the establishment of leading-edge vortex sheets and flow separation on such wings at 
high angles of attack; and (2) an increment in the aerodynamic moments, produced by the 
separated flow, which lags the motion of the configuration. The 0 derivatives are found 
to be very dependent on the frequency of oscillation, with the larger values obtained for 
the lower frequency. 

The paper also shows that the conventional use of rotary forced-oscillation data in 
the equation of motion to represent derivatives due to pure angular rates is erroneous at 
high angles of attack, where the 0 derivatives are of significant magnitude. In addition, 
it is shown that the conventional oscillating-airfoil theory and the lag-of-the-sidewash 
theory are inadequate for predicting the contribution of the vertical tail to Cn^. However, 
a flow-field-lag theory, devised in NACA RM L55H05 and extended herein, is found to yield 
values df C n ^ and which appear to be in qualitative agreement with experimental 

data for a current twin- jet fighter configuration. 

INTRODUCTION 

Recently, concern has arisen over the poor lateral-directional stability and control 
characteristics exhibited by many current high-performance military airplanes at high 

^Graduate research assistant, The George Washington University, Joint Institute for 
Acoustics and Flight Sciences. 



angles of attack. These poor characteristics have resulted in an alarming number of inad- 
vertent stalls, post-stall gyrations, and spins, which have produced a significant number 
of losses of aircraft and aircrews as well as severe operational restrictions. In view of 
this problem, the National Aeronautics and Space Administration is currently engaged in a 
broad research program designed to provide fundamental information regarding the aero- 
dynamic characteristics of fighter configurations at high angles of attack. It is intended 
that such information will serve as a basis for the design of future aircraft which will be 
either inherently spin resistant or which will use avionics and control- system concepts 
for automatic spin prevention. 

Research conducted by NACA in the 1950's (refs. 1 to 5) indicated that at high angles 
of attack, the magnitudes of the aerodynamic stability derivatives due to rate of change of 
sideslip p became quite large for swept- and delta-wing configurations, and that such 
derivatives had large effects on the dynamic stability of these configurations. Unfortu- 
nately, this research was conducted during a period of time when it was generally believed 
that future fighter aircraft would be standoff missile launchers, with no requirements for 
severe maneuvering or the attendant requirement for flight at high angles of attack. 

Because of an apparent lack of application of this research, the results of the past inves- 
tigations were fragmented and are therefore generally not considered in current analysis 
techniques for flight at high angles of attack. Recent military experience, however, has 
indicated a renewal of close-in, air-to-air combat involving strenuous maneuvering at 
angles of attack near the stall. The current stall/spin problems reflect this change in 
operational requirements, and it appears that a reconsideration of past research is 
required if future stall/spin problems are to be avoided. 

The purpose of the present report is to summarize, in some detail, pertinent infor- 
mation produced by NACA and NASA with regard to the derivatives due to rate of change 
of sideslip, and to emphasize the significance of these derivatives on the dynamic lateral- 
directional stability characteristics of high-performance swept-wing aircraft at high angles 
of attack. The report includes (1) a description of wind-tunnel test techniques used to 
measure the p derivatives; (2) illustrations of typical results; (3) a discussion of the 
effects of p derivatives on dynamic stability characteristics; and (4) a brief description 
of concepts which might be used to predict the values of such derivatives. 

A chronological listing of publications, not included in the reference section, which 
pertain to the low-speed, lateral-directional dynamic stability derivatives of aircraft is 
presented. 


SYMBOLS 

All longitudinal aerodynamic data are presented with respect to the wind system of 
axes. The lateral-directional aerodynamic data are presented with respect to the body 
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system of axes, as shown in figure 1, unless otherwise noted. Dimensional values are 
presented herein in the International System of Units (SI) with equivalent values given 
parenthetically in the U.S. Customary Units. 


a 


^m 

Cn 

C Y 

f d 


l xz 


nondimensional tail length, v - + 0.5| 

\ c v/2 ) 

wing span, m (ft) 

mean aerodynamic chord, m (ft) 

mean aerodynamic chord of vertical tail, m (ft) 

drag coefficient, Fj-j/q^S 

lift coefficient, Fjy/q^S 

rolling-moment coefficient, M^q^Sb 

pitching-moment coefficient, My/q^Sc 

yawing-moment coefficient, M2y / q oo Sb 

side-force coefficient, Fy^q^S 

drag force, N (lbf) 

lift force, N (lbf) 

side force, N (lbf) 

unsteady circulation function 

moment of inertia about longitudinal body axis, kg-m^ (slug-ft^) 
moment of inertia about normal body axis, kg-m^ (slug-ft^) 
product of inertia, kg-m^ (slug-ft^) 
unsteady circulation function 
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k reduced frequency of oscillation, cob/2V 

k v vertical-tail reduced frequency, a)C v y/2V 

l y distance from origin of axis to Cy/4, m (ft) 

M x rolling moment, m-N (ft-lbf) 

My pitching moment, m-N (ft-lbf) 

Mg yawing moment, m-N (ft-lbf) 

p rolling velocity, rad/sec 

q^ free-stream dynamic pressure, N/m^ (lbf/ft^) 

r yawing velocity, rad/sec 

S wing area, m^ (ft^) 

S y vertical-tail area, m^ (ft^) 

tj ^2 time required for oscillation to reach half-amplitude, sec 

V velocity, m/sec (ft/sec) 

y 0 amplitude of lateral oscillation, m (ft) 

a angle of attack, deg or rad 

j8 angle of sideslip, deg or rad 

j§ rate of change of sideslip angle, rad/sec 

X taper ratio 

ct sidewash angle (positive when effective sideslip at the vertical tail is 

reduced), deg or rad 

4>l,<p n phase angles associated with separation effects, rad 
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<j) frequency of oscillation, rad/sec 


( a Si 

increment of 

% 

due to vertical tail 

( 4 %) v 

increment of 

C Y 

Y 0 

due to vertical tail 

Subscripts: 





k quantity measured under oscillatory conditions 

s stability-axes data 

exp experimental values obtained under static conditions 

theo theoretical 
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BACKGROUND 


In the analysis of dynamic stability and control characteristics of airplanes, the 
lateral-directional stability derivatives due to j3 have traditionally been neglected because 
(1) the magnitudes of such derivatives were estimated to be small for the low values of 
angle of attack associated with the conventional flight envelope; and (2) little information 
was available regarding the magnitudes of the derivatives at high angles of attack. As will 
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be discussed, it was recognized that a lag of the sidewash at the vertical-tail location 
tended to increase aerodynamic damping in a manner similar to the lag of downwash for 
the longitudinal case; however, the lag of the sidewash did not appear to have large effects 
at low angles of attack for configurations under consideration at that time. 

As pointed out in reference 6, in the early 1950’s, several cases arose in which total 
disagreement existed between theoretical predictions of lateral-directional stability char- 
acteristics of certain swept -wing configurations and experimental results obtained during 
flight tests of dynamically scaled free-flight models and actual airplanes at high angles of 
attack. In these cases, theoretical methods generally predicted very unstable Dutch roll 
motions whereas flight tests showed highly damped stable motions. Subsequently, wind- 
tunnel test techniques devised to measure dynamic stability derivatives identified large 
values of 0 derivatives at high angles of attack, which, when used in the equations of 
motion, produced good agreement with flight tests. 

These wind-tunnel studies identified two principal sources of the 0 derivatives at 
high angles of attack: lag of sidewash at the vertical-tail location and lag of flow separa- 
tion and attachment on highly swept wings. In both situations, the time lags involved in 
flow over the vehicle resulted in large values of rolling and yawing moments which were 
not in phase with the motion. The derivatives became quite large at angles of attack near 
the stall, where sensitivity of the flow to small motions became severe and large areas of 
stalled flow were in evidence. 

WIND-TUNNEL MEASUREMENT TECHNIQUES FOR 0 DERIVATIVES 

Accurate analysis of the dynamic stability and control characteristics of airplanes 
depends directly on accurate aerodynamic inputs for computing time histories of aircraft 
motions and the stability of these motions. Theoretical prediction methods for the dynamic 
derivatives at high angles of attack have not advanced to the level of sophistication cur- 
rently reflected in such methods for low angles of attack. Most of the information per- 
taining to these dynamic derivatives has therefore been generated by specialized wind- 
tunnel test techniques. 


Rotary Forced-Oscillation Technique 

The most common wind-tunnel method in current use for obtaining dynamic stability 
derivatives is the rotary forced-oscillation technique illustrated in figure 2. In this tech- 
nique, a wind-tunnel model is forced to oscillate in roll or yaw at predetermined constant 
values of frequency and amplitude. A detailed discussion of the test setup and instrumen- 
tation required for a typical forced-oscillation technique is given in references 7 and 8. 
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It should be noted that because of kinematic constraints produced by this type of test, the 
measured dynamic stability parameters are a combination of the dynamic stability deriva- 
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tives. For example, forced-oscillation tests in yaw produce a parameter 
in body axes, or ^C nr j - ^C n ^j in stability axes. Unfortunately, the combined deriva^ 

tives cannot be separated or individually identified by using rotary forced- oscillation 
techniques. 


Linear Forced- Oscillation Techniques 

In addition to the rotary-type forced-oscillation technique, several devices have 
been employed to provide a linear sidewise oscillatory motion for the measurement of 
the /3 derivatives. The linear-oscillation technique permits the measurement of "pure" 

/§ derivatives over a wide range of angle of attack. A sketch of the setup involved for 
such tests is given in figure 3, and detailed descriptions of the apparatus and instrumenta- 
tion used may be found in references 2 and 4. 

Curved- and Rolling- Flow Techniques 

A unique wind tunnel was designed and constructed by the NACA at Langley in the 
early 1940’s to permit the measurement of derivatives due to "pure" angular rates. This 
tunnel, known as the Langley stability tunnel, was acquired by the Virginia Polytechnic 
Institute in 1958, and is currently in operation at that institute. The wind tunnel functions 
in a conventional manner for the measurement of static stability and control derivatives; 
however, it also incorporates interchangeable circular and square test sections which per- 
mit specialized tests for the measurement of dynamic derivatives. The circular test sec- 
tion is equipped with a motor- driven rotor which imparts a helical motion to the air stream 
and thus subjects a stationary model to a flow field similar to that which exists about an 
airplane in rolling flight. In this manner it is possible to measure aerodynamic deriva- 
tives due to "pure" steady-state roll rates. A more detailed discussion of the rolling-flow 
technique is given in reference 9. 

The square test section of the wind tunnel has a unique capability in that the vertical 
sidewalls are designed with sufficient flexibility so that they may be deflected into a curve, 
thus creating a curved airflow past the model. Jackscrews are positioned at regular inter- 
vals along each wall to allow the curvature to be set at predetermined values. In order to 
simulate flight in a curved path, it is necessary to redistribute the velocity profile in the 
radial direction in the tunnel. This need is accomplished by installing vertical wire 
screens in the flow upstream of the test section. These screens vary in mesh across the 
wind tunnel, with the most dense portion of the screens located nearer the center of curva- 
ture. A sketch showing a typical curved- flow test arrangement is shown in figure 4 and a 
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complete description of the tunnel and its operation is given in reference 10. With this 
test setup, it is possible to measure aerodynamic derivatives due to "pure" steady-state 
yaw rates. 


Pure-Yawing-Oscillation Technique 

In addition to the abovementioned techniques, an apparatus which superimposed both 
a yawing motion and a sidewise oscillatory motion was developed for the stability tunnel. 
With this apparatus, the resultant motion of the model was a pure yawing oscillation. A 
more detailed description of this test setup may be found in references 11 and 12. 

By using results obtained from the pure-yawing, rolling-flow, and rotary forced- 
oscillation tests, it is theoretically possible to add or subtract the magnitudes of the vari- 
ous test results and obtain values of the 0 derivatives. 

ILLUSTRATION OF TYPICAL RESULTS 
Isolated Wings 

Several investigations have been conducted to determine the physical causes of the 
j3 derivatives. In particular, tests conducted on isolated wings have provided some insight 
in this area. 

Static considerations .- Shown in figure 5 are an unswept wing, a 45° swept wing, and 
a 60° delta wing, which were investigated in reference 3. The static longitudinal charac- 
teristics of the wings are presented in figure 6(a). These data show that the unswept wing 
stalled at an angle of attack of about 16°, the swept wing stalled at an angle of attack of 
about 25°, but the delta wing was not stalled at the maximum angle of attack (30°) reached 
in the tests. The effective dihedral derivative C ^ as determined over the sideslip- angle 
range of ±10° is presented in figure 6(b) for each of the wings. The data show that the 
unswept wing experienced a large increase in positive dihedral effect near the stall, where- 
as the swept and delta wings exhibited large unstable values of (plpj near the stall. 

These characteristics are typical results for unswept and swept wings and may be inter- 
preted in a physical sense from the flow patterns indicated in figure 6(c). In addition, the 
mechanism which produces large 0 derivatives for isolated wings may be anticipated 
from these flow patterns. 

As shown in figure 6(c), the swept and delta wings shed strong vortices at high angles 
of attack, whereas the unswept wing does not. When subjected to sideslip, the advancing 
wing panel for the swept and delta wings becomes stalled, and there results a complete loss 
of dihedral effect. For the unswept wing the leading panel does not experience stall, but 
rather the trailing panel stalls, and a large increase in dihedral effect occurs as a result. 


8 



The foregoing characteristics are indicators of flow conditions which give rise to large 
values of the 0 derivatives during dynamic motions. 

Results of dynamic tests .- The basic phenomenon responsible for large 0 deriva- 
tives for swept and delta wings is the time lag required for the previously discussed vortex 
flows of these wings to follow the dynamic motions of the wing. This characteristic has 
been shown for a delta wing in reference 13. In a static condition at high angles of attack 
and zero sideslip, the delta wing sheds two strong vortices which emanate from the apex of 
the wing and remain symmetrically displaced over the wing. Motion-picture photographs 
presented in reference 13 show that when the wing is subjected to oscillatory motions, the 
vortex flow field lags behind the model motion. Because of the flow lag the vortices are 
asymmetrically displaced at a time when the model is at zero sideslip. The resulting flow 
field thus creates a rolling moment which is due to the rate of change of sideslip. The 
foregoing discussion is a brief description of the mechanism wherein large values of 
are produced by swept and delta wings at high angles of attack. 

Examples of the 0 derivatives measured for a 60° delta wing and a 45° swept wing 
are reported in reference 5 and presented in figure 7. The static longitudinal characteris- 
tics of the wings (tested in refs. 5, 11, 12, and 14) are presented in figure 8 and the static 
lateral-directional characteristics are presented in figure 9 for reference. Presented in 
figures 10 and 11 are representative data pertaining to the 0 derivatives for these wings. 

Figure 10 gives a comparison of data for the delta wing obtained by the forced- 
oscillation-in-yaw technique (ref. 14) and the linear -oscillation technique (ref. 5) at two 
values of reduced frequency k. Also shown are data obtained from steady-state yawing 
tests (ref. 11) and the pure-yawing-oscillation tests (ref. 12). It should be noted that the 
results from the steady-state yawing tests and the pure-yawing-oscillation tests are in 
disagreement. The reason for the discrepancy is unknown; however, it should be noted 
that for either case the magnitudes of the 0 derivatives are larger than those of the 
derivatives due to pure yawing at high angles of attack, and the results of the forced- 
oscillation test technique provide a qualitative measure of these derivatives for isolated 
wings. The magnitudes of the derivatives are seen to be highly dependent on the value 
of k, with the larger values obtained for the lower value of k. 

Similar results obtained for the 45° swept wing are presented in figure 11. The data 
show trends similar to those exhibited by the delta wing; however, agreement between 
results of the forced-oscillation and linear-oscillation techniques is not as good as the 
agreement obtained for the delta wing. In both cases it should be noted that the results 
indicate large values of and (Clfcj 5 311(1 1<; ma Y be anticipated that in equations 

of motion, the use of forced-oscillation results as "pure" yawing derivatives would be com- 
pletely erroneous for the configurations tested. 
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It might be expected that the magnitude of the /3 derivatives near the stall would be 
highly dependent on the amplitude and frequency of motion used in the forced-oscillation 
test technique. Although not presented herein, the results of several studies indicate that 
these derivatives may have large dependence on amplitude and frequency, and large non- 
linear effects may be present at high angles of attack. In addition, because of the nature 
of the flow field associated with the 0 derivatives, Reynolds number effects may be 
significant. 

The past research has also indicated that the 0 derivatives of unswept wings may 
be much smaller than those of swept wings at high angles of attack. For example, forced- 
oscillation test results (from ref. 3) for the unswept, the 45° swept, and the 60° delta wings 
of figure 5 are presented in figure 12 for two values of k. These data show that the mag- 
nitudes of |Cn r ^ - and (^Irj ~ are muc h smaller for the unswept wing 

s s 

at high angles of attack. This result serves to emphasize the importance of wing sweep 
in the formation of 0 derivatives at high values of angle of attack. 

Contribution of Vertical Tail 

The second major contributor to the 0 derivatives at high angles of attack is the 
vertical tail. This contribution is a result of the lag of sidewash at the vertical-tail loca- 
tion. The lag of the sidewash, and its effect on the damping in yaw of the tail, is analogous 
to the lag of downwash and its effect on the damping in pitch of a horizontal tail, which was 
first discussed by Cowley and Glauert in reference 15. The lag of downwash was treated 
in the past as an additional angle of attack of the horizontal tail which was due to the time 
required for the wing disturbance to travel the distance between the wing and the horizon- 
tal tail. The case for the lag of sidewash at the vertical tail follows from the same phys- 
ical mechanism, but the characteristic length involved is not as well defined. 

Typical results of forced-oscillation-in-yaw tests for fighter configurations at high 
angles of attack indicate that the vertical tail may add considerable damping at high angles 
of attack, even though the tail may make no significant contribution to static directional 
stability. For example, shown in figure 13(a) are results obtained during static wind- 
tunnel tests for a current fighter configuration which is subject to loss of directional 
stability at high angles of attack. The data show that the contribution of the vertical tail 
to C n ^ was markedly reduced above a = 30°. This result is usually attributed to a 
combination of reduced dynamic pressure and adverse sidewash at the vertical-tail loca- 
tion. However, as shown in figure 13(b), the combination of derivatives measured during 
forced-oscillation-in-yaw tests shows an increased vertical- tail contribution at high angles 
of attack. As indicated previously, this result may be attributed to the vertical-tail con- 
tribution to C n ^, which is associated with a lag of the sidewash. Analytical and exper- 
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imental results for the vertical-tail contribution to C n * are discussed in a subsequent 

• P 

section '^Prediction Methods for |3 Derivatives." 

EFFECT OF /3 DERIVATIVES ON DYNAMIC LATERAL-DIRECTIONAL 
STABILITY CHARACTERISTICS 

The results of reference 1 show that the inclusion or omission of the fi derivatives 
in theoretical equations of motion can produce large differences in the calculated dynamic 
stability characteristics. In particular, reference 1 shows that when data obtained from 
forced-oscillation tests — which represent a combination of derivatives — are used as a 
total value for the pure angular-rate derivatives and the jS derivatives are assumed to 
equal zero, the results vary significantly from the results obtained when all pure deriva- 
tives are used. In addition, reference 16 has shown that the frequency effects of stability 
derivatives can cause considerable changes in predicted aircraft motion. In order to 
illustrate these results, calculations are made for a hypothetical delta-wing fighter 
configuration. 

The mass and inertial properties used in the calculations are presented in table I and 
the aerodynamic characteristics for the configuration, obtained from references 4, 17, and 
18, are presented in figure 14. It should be noted that the static lateral- directional char- 
acteristics of this configuration (fig. 14(b)) are typical of those exhibited by many current 
fighter configurations. In particular, the configuration exhibits a marked reduction in both 
directional stability and dihedral effect at high angles of attack. Experience has shown 
that these conditions usually result in a directional divergence, or "nose slice." The 
yawing and rolling derivatives (fig. 14(c)) were obtained by the rolling- and curved-flow 
test techniques previously discussed. The derivatives due to rate of change of sideslip 
(fig. 14(d)) were obtained from wind-tunnel tests in reference 4, and the magnitude of Cy^ 
was assumed to be zero for the calculations. The dynamic lateral- directional stability 
characteristics were calculated by means of classical three-degree-of-freedom linearized 
equations for the configuration in trimmed flight at angles of attack of 20° and 28° at an 
assumed altitude of 7620 m (2 5 000 ft). The condition of a = 20° represents a case for 
which the static and dynamic data indicate a stable configuration with relatively small 
values of the j3 derivatives. At a = 28°, the configuration exhibits static instability 
(negative value of Cn^> a marked reduction in magnitude of C ^ , and large values of 
the /3 derivatives. 

The results of the calculations are presented in table II in terms of the time to halve 
amplitude tj/2 and the nondimensional reduced frequency k of the various lateral- 
directional modes of motion. Positive values of t represent damped (dynamically 
stable) modes of motion, whereas negative values represent undamped (dynamically 
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unstable) modes of motion. Case 1 represents the results of calculations in which the 
true values of the angular-rate derivatives were used and the 0 derivatives were omit- 
ted. Cases 2 to 4 represent the results of calculations in which the angular-rate deriva- 
tives were used, and the experimental values of the 0 derivatives based on data obtained 
for three different values of k were included. Cases 5 to 7 represent results for which 
the 0 derivatives were summed with the pure rate derivatives to form combinations 
similar to those obtained from conventional forced-oscillation tests. For these cases, the 

9 

resulting sums were used as values for the rate derivatives and the 0 terms in the 
equations were set equal to zero. 

In table n(a), the results obtained from the calculations for a = 20° are presented. 
The results show that including the 0 derivatives (cases 2 to 4) increased the damping of 
the Dutch roll mode but had essentially no effect on the frequency of the Dutch roll mode 
or the damping of the roll and spiral modes. Combining the derivatives (cases 5 to 7) had 
little additional effect on the Dutch roll and roll modes; however, the spiral mode became 
unstable. The foregoing results, obtained for relatively small values of the 0 deriva- 
tives, indicate relatively minor effects caused by neglect or misuse of the 0 derivatives 
at low angles of attack. For a = 28° (table n(b)), the data show pronounced differences 
caused by the relatively large 0 derivatives. Use of the angular-rate derivatives and 
omission of the 0 derivatives (case 1) resulted in a stable Dutch roll mode, a virtually 
neutral stable roll mode, and a very unstable spiral mode. However, when the 0 deriva- 
tives were included in the calculations (cases 2 to 4) the Dutch roll mode ceased to exist, 
and two additional aperiodic modes were formed. In addition, it is seen that the magnitude 
of the 0 derivatives used in the calculations has an effect on the computed values of the 
time to halve amplitude of the spiral and aperiodic modes. When the angular-rate and 0 
derivatives were used in combination (cases 5 to 7), the damping of the entire system was 
redistributed, resulting in highly stable spiral and roll modes and a highly unstable Dutch 
roll mode. The foregoing results illustrate that in equations of motion, the conventional 
use of the rotary forced-oscillation data to represent derivatives due to pure angular 
rates is erroneous at high angles of attack, where the 0 derivatives are of significant 
magnitude. 
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In addition to having effects on dynamic stability, it would be expected that the 0 
derivatives would also have a large effect on the design of control systems for flight con- 
ditions at high angles of attack and on related analysis techniques. For example, a recent 
study (ref. 19) has shown that omission of the 0 derivatives may produce considerable 
errors in parameter-identification techniques. 

PREDICTION METHODS FOR 0 DERIVATIVES 

The previous discussion has emphasized the large effects caused by conventional use 
of forced-oscillation data in equations of motion at high angles of attack. Unfortunately, 
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no analytical estimation method currently exists for prediction of the /3 derivatives, and 
the subject has received little attention in the literature. At the very least, a method is 
required to permit a realistic qualitative separation of the combined derivatives measured 
in conventional forced-oscillation tests. Such methods should consider the contributions 
of both the vertical tail and the wing to the j3 derivatives. 


Vertical-Tail Contribution 

Oscillating-airfoil theory .- Expressions for the lift and pitching moment arising 
from the unsteady circulation about a two-dimensional airfoil subjected to small sinusoidal 
oscillations were developed in reference 20. This analysis was extended to include the 
case of finite- span wings in reference 21 and applied to the case of an oscillating vertical 
tail in references 22 and 23. From the methods presented in references 22 and 23, the 
following expression for ^AC n ^ is developed: 




Sy 

S 


+ a + 



(1) 


where the unsteady circulation functions J and G are tabulated in reference 23 as a 
function of the aspect ratio and reduced-frequency parameter k v of the vertical tail. 
Although the oscillating-airfoil theory is frequency dependent, its application to the present 
subject is extremely limited because it does not consider sidewash or separated-flow 
effects. 

Lag of sidewash .- The following expression for ^AC n ^j is developed in refer- 
ence 24 and is based on the hypothesis that the sidewash at the vertical tail affects the 
yawing moment through a time lag t, where r * Z v yV: 



(The negative sign has been introduced for consistency with the present notation for side- 
wash.) Essentially, this development assumed that the sidewash originates at the center 
of gravity of the configuration and requires a finite time for the flow field to impinge on 
the vertical tail. For a simple fuselage — vertical-tail model, with auxiliary fins placed at 
the center of gravity to generate sidewash, reference 24 shows that equation (2) yields 
results which correlate fairly well with experimental data. 

Figure 15 shows the vertical-tail contribution to C n ^ for the delta-wing fighter 
configuration of reference 4. The vertical-tail contribution was determined by subtracting 
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the experimental wing-fuselage data from the data for the complete configuration. Fig- 
ure 15 also presents the vertical-tail contribution calculated with the oscillating- airfoil 
theory (eq. (1)) and the lag-of-the-sidewash theory (eq. (2)). In these calculations, the 
geometric characteristics were obtained from table I, and the unsteady circulation func- 
tions were obtained from reference 23; the increment (aCy^J due to the vertical tail 

and the sidewash parameter 8cr/8j3 were obtained from reference 17. 

From figure 15 it is seen that both the theories are completely inadequate for pre- 
dicting the vertical-tail contribution to C n ^. The lack of correlation of the oscillating- 
airfoil theory is probably due to the neglect of the sidewash generated by the wings; the 
lack of correlation of the lag-of-the-sidewash theory is probably due to the fact that the 
sidewash at the vertical tail is not necessarily generated at the center of gravity, as was 
assumed in the development of that theory. 

Wing Contribution 

As discussed previously, reference 13 has shown that under oscillatory conditions 
the flow field of a 60° delta wing lags behind the model motions. This lag has been inves- 
tigated in detail (see ref. 3) by means of oscillograph traces of the output of the yawing- 
moment channel of a strain-gage balance. Additionally, reference 3 considers the discrep- 
ancy between the theoretical static lateral- directional stability derivatives C n ^ and 
(see ref. 25) and those measured in both conventional static force tests and in forced- 
oscillation tests. On the basis of these considerations, reference 3 has shown that the 
large values of ^C nr j - (CnpJ measured for delta wings at high angles of attack can be 

attributed to an incremental yawing moment which is produced by separated flow and which 
lags the motion by a roughly constant time interval. In addition, the increase in the mag- 
nitude of - (CnpJ increasing angle of attack is not caused by an increase 

in the time lag and is therefore probably attributable solely to the increased strength of the 
vortex flow and separation at the higher angles of attack. This result has led to the devel- 
opment of the following expression for C n ^ ( see re f. 3); 





sin <t> n 


exp 


(3) 


where 0 n is the phase angle associated with separation effects. 

Since equation (3) is presented in reference 3 without a detailed development, it is 
derived in the appendix. A similar expression for • is obtained: 
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( 4 ) 



where <p^ is the phase angle associated with separation effects. Additionally, since the 
values of cp n and are obtained by examination of oscillograph traces in reference 3, 
expressions for and are developed in the appendix and presented below: 


and 



(5) 


( 6 ) 


The use of the flow-field-lag theory has been shown in reference 3 to yield results 
which correlate well with experimental results of a 60° delta wing over the range of fre- 
quency considered. 


As an illustration of the applicability of this technique to complete configurations, 
the twin-jet fighter airplane studied in references 7, 26, and 27 is considered. Figure 16 
shows the variation of the lateral-directional stability derivatives C n ^ and with 

angle of attack. These data were obtained in reference 7 from conventional static force 
tests and from the in-phase data obtained for forced-oscillation-in-yaw tests ^assuming 
that k C n * and k are negligible) at a reduced frequency k of 0.156. Also pre- 
sented in figure 16 are approximate theoretical values of C n ^ and C j a . These approxi- 


'h S' 


mate values are developed by first obtaining the theoretical variation of and 

(Ci^j with lift coefficient for the wing-fuselage combination. (See ref. 25.) These values 

were then transferred to the body-axis system, and the experimentally determined incre- 
ments to Cnp and due to the tail (evaluated at a = 0) were added. The approxi- 
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mate theoretical values therefore represent reference values which would be obtained if 
separation effects were not present. 


Figure 17 presents the calculated values of the acceleration-in- sideslip derivatives 
C n ^ and C 'i* obtained by applying the flow- field-lag theory (eqs. (3) to (6)) to the data of 
figure 16. Also presented in figure 17 are values of Cn^ and obtained by subtrac- 

tion of the data measured under conditions of pure yawing flow (ref. 27) from the data 
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measured with the forced-oscillation-in-yaw technique (ref. 7) at a reduced frequency of 
oscillation k of 0.156. Figure 17 shows that the flow-field-lag theory yields values of 
C n ^ and which are in reasonable agreement with the data obtained from the com- 

bination of tests with yawing flow and forced oscillation in yaw. Although the foregoing 
method could be carried out for the present configuration only over a limited angle-of- 
attack range because equations (5) and (6) yield values of the cosine function greater than 
unity, it offers some promise for future development of predictive techniques. 

CONCLUSIONS 


From the present summary of experimental and theoretical results pertaining to the 
aerodynamic stability derivatives due to the rate of change of sideslip /3 the following 
conclusions are made: 

1. The 0 derivatives and C^ are large for swept- and delta-wing con- 

figurations at high angles of attack. 

2. The physical flow phenomenon responsible for the 0 derivatives is associated 
with the establishment of leading- edge vortex sheets and flow separation on the wings at 
high angles of attack. Additionally, the large values of the 0 derivatives can be attrib- 
uted to an increment in the aerodynamic moments produced by the separated flow, which 
lags the motion of the configuration. 

3. The /3 derivatives are very dependent on the frequency of oscillation, with the 
larger values obtained for the lower frequency. 

4. In equations of motion, the conventional use of the rotary forced-oscillation data 
to represent derivatives due to pure angular rates is erroneous at high angles of attack, 
where the 0 derivatives are of significant magnitude. 


5. Both the oscillating-airfoil theory and the lag-of-the sidewash theory are inade- 
quate for predicting the contribution of the vertical tail to C n ^. 

6. The flow- field- lag theory, devised in NACA RM L55H05 and extended herein, is 
' found to yield results which are in reasonable agreement with experimental data for a 

current twin- jet fighter aircraft. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
June 17, 1975 
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APPENDIX 


DEVELOPMENT OF EQUATIONS FOR C n ^ AND 
BASED ON FLOW- FIELD- LAG THEORY (REF. 3) 


The equations developed herein for C n ^ and C ^ are based on the flow- field- lag 
theory (ref. 3) and are developed for the case of a pure sidewise oscillation. However, 
as will be discussed, they are applicable to the data obtained under conditions of forced 
oscillation in yaw. 

For the pure- sidewise- oscillation technique, discussed in references 2 and 4, the 
yawing velocity and yawing acceleration are identically zero and the model undergoes 
continuous changes in sideslip and rate of change of sideslip. Therefore, the yawing- 
moment signal C n (t) may be expressed as 


c n ,t> . + (c» .)kt) 

where, for harmonic oscillation (see ref. 4), 
y w 

/3(t) = cos a)t 

Introducing the nondimensional frequency parameter k = u>b/2V yields 
2y 

/3(t) = -r-2 k cos cot 
b 


(Al) 


(A2) 


(A3) 


Differentiating equation (A3) with respect to time and introducing the nondimensional 
form of /3 (that is, /§ = /3b/2V) yields 


• 2y^ 0 

j8 = - — ^ k^ sin cut 
b 


(A4) 


Substituting equations (A3) and (A4) into equation (Al) yields 

c n(t) = ( C n^) k "E 2 k cos wt " ( Cn ^) k IT k2 Sin 

Multiplying equation (A5) by cos u>t and integrating over the period yields 

. p2ir 

(S) k = 2^J 0 C n« COSwt " dt 


(A5) 


(A6) 
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Similarly, multiplying equation (A5) by sin cot and integrating over the period yields 

b C 2it 


(%) 


k 2y k it 


(*2ir 

JO Cn(t) 


sin cot co dt 


(A7) 


Adding and subtracting a theoretical value of C n ^ (which would be obtained in the 
absence of flow-separation effects, see ref. 25) to equation (A6) yields 


(sl ■ N theo + £ c - (t) cos - (s) 


(A8) 


theo 


Since the (3 derivatives are assumed to arise from flow- separation effects, a 
theoretical yawing-moment coefficient may be developed in terms of the theoretical values 
c np‘ This theoretical yawing-moment coefficient assumes no flow separation and may 
be written as 


M th e 0 SSL„^ C08ut 

Multiplying equation (A9) by cos cot and integrating over the period yields 

fc n \ * Tr~i — f 27r |Cn(t)l u cos cot co dt 

V p ) theo 2y o k7rJ0 L nV ^theo 

Substituting equation (A10) into equation (A8) yields 

b r 2ir 


(A9) 


(A10) 




l + So ' { c » (t) - C c - (t 3 theo } - -* - * < A11 > 


Assuming 


:,l(t) ' ^theo ‘ KU ‘ (S) 


i 2 y 0 k 


theo, 


cos 


(cot - <£ n ) 


(A12) 


where 0 is the phase angle associated with separation effects, and substituting equa- 
tion (A12) into (All) yield, upon carrying out the integration, 


w. 


= (c- 


-n. 


theo 


Nexp ' N 


theo 


cos 0 n 


(A13) 
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Rearranging equation (A13) and solving for <fi n yields 

(SLo ' (%) k 


*n = 


COS 


(Stheo'lS) 


exp 


(A14) 


Substituting equation (A9) into equation (A12) and substituting the resulting expression for 
C n (t) into equation (A7) yield upon carrying out the integration, 


Cn /3 = 


Ntheo ' N«p 


sin <f) n 


(A15) 


where <£ n is given by equation (A14). 

In a similar fashion it can be shown that 


(%L '(%) 

' 7 theo ' H/ exp 


sin (p 7 


(A16) 


where is the phase angle associated with separation effects and is given by the 
expression 




= cos 


-1 


(%), 


theo 


% 


( C ^theo ’ ( l ») 


exp 


(A17) 


In the development of these equations, a sideward oscillatory motion was assumed; 
however, it should be noted that the expressions for C n £ and are functions of the 

in-phase derivatives and . Although the in-phase data obtained from the 

forced-oscillation-in-yaw tests are a combination of derivatives, such as C nfi + k^C n . 

2 P ^ 

and + k C^j,, experience has shown that the magnitudes of the r derivatives are 

usually negligible and the total in-phase components may be considered as pure derivatives 
(C n/3 ) k and (%) k - Therefore, the method in this appendix for evaluating C n £ and 

is also applicable to the data obtained from the forced-oscillation-in-yaw technique. 
Additionally, it should be noted that there is no restriction as to the use of the stability or 
body system of axes. 
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TABLE L- MASS AND GEOMETRIC CHARACTERISTICS OF THE 
DELTA- WING AIRPLANE 


Mass, kg (slugs) 

Moments of inertia, kg-m 2 (slug-ft^): 



X XZ 

Overall length, m (ft) 

Wing: 


Span, m (ft) 

Area, m 2 (ft 2 ) 

Mean aerodynamic chord, m (ft) 

Aspect ratio 

Taper ratio 

Dihedral, deg 

Incidence, deg 

Leading-edge sweep, deg .... 
Airfoil section 

Vertical tail: 

Span, m (ft) 

Area, m 2 (ft 2 ) 

Root chord, m (ft) 

Taper ratio 

Leading- edge sweep, deg 

Airfoil section 

Mean aerodynamic chord, m (ft) 

Area ratio, Sy/S 

Tail-length ratio, l v f b 

Aspect ratio 

Tail length, m (ft) 


. 17 727 (1215) 

40 606 (29 950) 
229 859 (169 538) 
7106 (5241) 

. 15.90 (52.16) 

. 10.75 (35.26) 

. 50.00 (538.30) 
6.21 (20.38) 

2.31 

0 

0 

0 

60 

. . NACA 65A003 

. . 3.53 (11.59) 

. . 5.72 (61.59) 

. . 3.24 (10.63) 

0 

42.5 

. . NACA 65-006 
. . . 2.16 (7.10) 

0.115 

0.59 

2.18 

. . 6.33 (20.77) 
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TABLE H.- CALCULATED DYNAMIC STABILITY CHARACTERISTICS 


(a) a - 20° 

Aerodynamic derivatives Results 



a Reduced frequency at which /3 derivatives were measured. 
b j3 derivatives combined with pure angular- rate derivatives. 


(b) Oi se 28° 

Aerodynamic derivatives Results 

I | v T 1 ' ' -r* 

Dutch roll Roll Spiral Aperiodic 


Case k mode mode mode modes 



a Reduced frequency at which j3 derivatives were measured. 

CO h • 

co (3 derivatives combined with pure angular-rate derivatives. 










00 

4 ^ 




Type Delta 

Sweep 60°(LE. ) 

Aspect ratio 2. 31 

X 0 

Airfoil NACA65 - 006.5 

Figure 5.- Sketch showing geometric characteristics of wings tested in reference 3. 
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K£KI Separated flow 
I I Attached flow 




Low angles of attack 


High angles of attack 


(c) Upper-surface flow characteristics in sideslip. 


Figure 6.- Concluded. 










(a) Low frequency. (b) High frequency. 

Figure 10.- Comparison of data obtained by various test techniques for 60° delta wing, 














a, deg 

(a) Static directional stability derivative. 



a, deg 

(b) Damping in yaw from forced-oscillation tests. 


Figure 13.- Effect of vertical tail on yawing derivatives of a current fighter configuration. 
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(d) Derivatives due to rate of change of sideslip (from ref. 4) 
Figure 14.- Concluded. 
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a, deg 

Figure 16.- Variation of static lateral-directional stability derivatives with angle 

of attack for twin-jet fighter airplane. 
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